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Distributions of event-by-event fluctuations of the mean transverse momentum and mean trans- 
verse energy near mid-rapidity have been measured in Au+Au collisions at ^/sjvjv = 130 GeV at 
RHIC. By comparing the distributions to what is expected for statistically independent particle 
emission, the magnitude of non-statistical fluctuations in mean transverse momentum is determined 
to be consistent with zero. Also, no significant non-random fluctuations in mean transverse en- 
ergy are observed. By constructing a fluctuation model with two event classes that preserve the 
mean and variance of the semi-inclusive pr or er spectra, we exclude a region of fluctuations in 
^/sjvjv = 130 GeV Au+Au collisions. 
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I. INTRODUCTION 

Phase instabilities near the QCD phase transition can 
result in non-statistical fluctuations that are detectable 
in final state observables Q]. These instabilities, which 
may occur due to random color fluctuations critical 
behavior at the QCD tri-critical point H, or fluctuations 
from the decay of a Polyakov loop condensate M, can 
result in a broadening of the transverse momentum or 
transverse energy distributions of produced particles for 
different classes of events. This phenomenon is expected 
to be detected experimentally by searching for deviations 
of the distributions of the event-by-event mean transverse 



momentum M PT or mean transverse energy M eT of pro- 
duced particles from the random distributions expected 
for statistically independent particle emission. 

An event-by-event analysis of M PT was previously per- 
formed for 158 A GeV/c Pb+Pb Collisions at the CERN 
SPS by Experiment NA49 [0. In that analysis, the M PT 
distributions measured over the rapidity range 4 < y v < 
5.5 and px range 0.005 < p? < 1.5 GeV/c, were found to 
be consistent with random fluctuations. NA49 also per- 
formed an event-by-event analysis of the K/n ratio ml, 
showing only very small deviations from random fluctu- 
ations. With an increase of ^/sjvjv to 130 GeV in RHIC 
collisions, unprecedented energy densities have been ob- 
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served [|10| , hence conditions may be more favorable for a 
phase transition from hadronic matter to a Quark-Gluon 
Plasma which may be indicated in non-random fluctu- 
ations. Presented here is an event-by-event analysis of 
M PT fluctuations and the first measurement of M eT fluc- 
tuations at mid-rapidity at RHIC. 

II. ANALYSIS 

The PHENIX experiment Q| consists of four spectrom- 
eters designed to measure simultaneously hadrons, lep- 
tons, and photons produced in nucleus-nucleus, proton- 
nucleus, and proton-proton collisions at RHIC. The two 
central arm spectrometers, which are located within a 
focusing magnetic field and each cover ±0.35 in pseudo- 
rapidity and A</> = 90° in azimuthal angle, are utilized 
in this analysis. The primary interaction trigger was de- 
fined using the Beam-Beam Counters (BBC) || and Zero 
Degree Calorimeters (ZDC) @. Events are selected with 
a requirement that the collision vertex along the beam 
axis has \z\ < 20 cm as measured by both the BBC and 
ZDC. Event centrality is defined using correlations in the 
BBC and ZDC analog response as described in [jl0|. For 
the present analysis, the events are classified according 
to the - 5%, - f 0%, 10 - 20%, and 20 - 30% most 
central events. 

The drift chamber [ pT| is used in conjunction with 
the innermost pad chamber, called PCI, to measure 
the transverse momentum of charged particles travers- 
ing the PHENIX acceptance. A fiducial section of the 
drift chamber is chosen to minimize the effect of time- 
dependent variations in the performance of the detector 
during the data-taking period. The fiducial volume of the 
M PT analysis spans an azimuthal range of A(j) = 58.5° 
and covers the pseudorapidity range \rj\ < 0.35. Recon- 
structed tracks [fl2f are required to contain a match to 
a hit in PCI to insure that the tracks are well recon- 
structed in three dimensions for reliable momentum de- 
termination. 

The M eT distribution is determined from clusters re- 
constructed in the two instrumented sectors of the Lead- 
Scintillator electromagnetic calorimeter [[7| Jl^ , p^| . The 
quantity ex is defined as the transverse energy per recon- 
structed calorimeter cluster as described in [O], which 
can include clusters that have been merged. The effects 
of cluster merging on the M eT distribution are discussed 
later. The fiducial volume of the M eT analysis spans an 
azimuthal range of A<fi = 45° and covers |ry| < 0.35. 

There are no acceptance nor efficiency corrections ap- 
plied to the semi-inclusive px or ex distributions prior 
to the calculation of M PT or M eT . These corrections do 
not vary from event to event and are identical for data 
and mixed events; therefore they do not modify the val- 
ues of the fluctuation quantities defined later. The M PT 
distributions are calculated using the formula 



Ntrachs 

M PT = (l/N tracks ) Pth (2- 1 ) 

i=i 

where N tra cks is the number of tracks in the event that 
pass the cuts outlined above and lie within the px range 
0.2 < pt < f.5 GeV/c. Similarly, the M eT distributions 
are calculated using the formula 

NduB 

M eT = (1/N clus ) e Tz, (2.2) 

i=l 

where N c i us is the number of calorimeter clusters in the 
event that lie within the ex range 0.225 < ex < 2.0 GeV. 
An event is excluded from the analysis if N tra cks or N c i us 
is below a minimum value to insure that there are a suf- 
ficient number of tracks or clusters to determine a mean 
and to exclude background events. This minimum value 
for the - 5%, - f 0%, f - 20%, and 20 - 30% central- 
ity classes, respectively, is 40, 30, 20, and 10 for the M PT 
analysis and 30, 20, f 0, and f for the M eT analysis. Ta- 
ble I lists statistics pertaining to the data samples used 
to determine M PT and Table H lists the statistics per- 
taining to the data samples used to determine M eT . The 
events used for the M PT and M eT analyses are considered 
independently of each other. 

In order to compare the M PT and M eT distributions 
to what is expected for statistically independent particle 
emission, mixed events are considered as the baseline for 
the random distribution. To obtain a precision compar- 
ison, it is important to match the number of tracks or 
clusters along with the mean of the semi-inclusive distri- 
bution of the mixed events to that of the data. There- 
fore, in both analyses, mixed events are constructed by 
pre-determining the number of charged particle tracks or 
calorimeter clusters in the mixed event N m i X by directly 
sampling the corresponding data N tra cks or N c i us distri- 
butions. Figure f shows a comparison of the N tra cks dis- 
tributions from the data and the normalized mixed event 
N m ix distribution for the — 10% centrality class. Once 
Nmix is determined, a mixed event is filled with px or 
ex values from the data with the following criteria: a) no 
two px or ex values from the same data event are allowed 
to reside in the same mixed event, b) only px or ex val- 
ues passing all cuts in the determination of M PT or M eT 
from the data events are placed in a mixed event, and c) 
only data events from the same centrality class are used 
to construct a mixed event corresponding to that class. 
Once a mixed event is filled with N m i X tracks or clusters, 
its M PT or M eT is calculated in the same manner as for 
the data events. 

For both analyses, the data contain a fraction of tracks 
or clusters within close physical proximity that have 
merged into a single track or cluster. This fraction is 
estimated by embedding simulated single-particle events 
that are processed through a detailed simulation of the 
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detector response into real data events, which are then 
reconstructed in the same manner as the data. For the 
— 5% centrality class, we estimate that 6% of the tracks 
and 5% of the clusters are affected. 

For the M PT analysis, tracks that are merged into a 
single reconstructed track typically have similar values of 
Pt- The result is a slightly lower value of N tra cks which 
causes a slight broadening in the width of the M PT distri- 
bution due to the reduced statistics per event. However, 
since the N tra cks data distribution is directly sampled 
during the construction of mixed events, the effect of 
merged tracks cancels for comparisons between the data 
and mixed events. 

For the M eT analysis, the effect of merged clusters is 
complicated by the fact that a single cluster is recon- 
structed with an er corresponding to the sum of the two 
(or more) particles contributing to the cluster. To under- 
stand this effect on the mixed events, we note that the 
fraction of merged clusters within a data event increases 
with event multiplicity. Also, many of the data events 
with the lowest M eT coincide with the lowest multiplic- 
ity events since they contain few, if any, merged clusters 
that would yield a higher M eT . When the merged clusters 
in the data events are randomly redistributed among the 
mixed events, low multiplicity mixed events can contain 
more merged clusters than the data events with the same 
multiplicity, resulting in a gross upward shift in M eT for 
those mixed events. This results in apparent excess non- 
random fluctuations at low M eT . Conversely, high mul- 
tiplicity mixed events can contain fewer merged clusters 
than the data events with the same multiplicity, result- 
ing in a gross downward shift in M eT for those mixed 
events. However, since the mean is taken over more clus- 
ters in this case, the effective shift in M &T is reduced at 
high M eT , and the apparent non-random fluctuations are 
much less pronounced. An estimate of the magnitude of 
this effect is presented later. 



III. RESULTS 



where 



(< p\ > ~ < pr > 
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The reciprocal of b is the inverse slope parameter of 
the pt distribution. With the track multiplicity distri- 
bution assumed to be a negative binomial distribution, 



fNBD(N t racks, 1/fc, < Ntracks >), 

can be calculated using 



the M PT distribution 
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9(1>t)= £ fNB D (N,l/k,<N >)f r (pT,Np,Nb), 

N=N min 

(3.4) 

where the loop is over N tra cks from N m i n to N max , which 
are the limits of the multiplicity. The value of the nega- 
tive binomial distribution parameter k is given by 



Pt 
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^ -^tracks ^ ^tracks 
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Therefore, given < pt >, <J PT , and < N tra cks > extracted 
from the semi- inclusive pt distribution, the expected ran- 
dom M PT distribution can be calculated. Figure 2 shows 
the M PT distribution for the 0—5% centrality class. Over- 
layed on the data as a dotted curve is the result of the 
calculation. The agreement between the data distribu- 
tion and the calculation illustrates the absence of large 
non-statistical fluctuations in the data. The remainder 
of this article will quantify the amount of non-statistical 
fluctuations observed and place limits on the level of fluc- 
tuations that can be present in central Au+Au collisions 
at y/sNN = 130 GeV. 

To quantify the magnitude of the deviation of fluctu- 
ations from the expectation of statistically independent 
particle emission, the magnitude of the fluctuation lot in 
the transverse quantity Mr, representing M l 
is defined as 



PT or M eT , 
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< M T > < M T > ' 
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To compare directly the semi-inclusive pt distribution 
to the M PT distribution assuming statistically indepen- 
dent particle emission, the closed form prescription out- 
lined in fl5|j is used. This prescription describes the semi- 
inclusive pt distribution using a Gamma distribution, 



The value of lot is calculated independently for the data 
distribution and for the baseline, or mixed event, distri- 
bution. The difference in the fluctuation from a random 
baseline distribution is defined as 



W (T, data) ~ W (T, baseline)- 



(3.7) 



f(pr) = fr(PT,P,b) 



(bp 7 



l e -bpi 



(3.1) 



where p and b are free parameters that are related to 
the mean and standard deviation of the semi-inclusive 
distribution as 



< Pt >' 



<Pt > 



(3.2) 



The sign of d is positive if the data distribution contains a 
correlation, such as Bose- Einstein correlations |ll|, when 
compared to the baseline distribution. The fraction of 
fluctuations which deviate from the expectation of sta- 
tistically independent particle emission is given by 
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where crrr, data) refers to the standard deviation of the 
event-by-event Mr data distribution and <j(r, baseline) is 
the corresponding quantity for the baseline, or mixed 
event, distribution. In the absence of a common lan- 
guage for the analysis of M PT and M £T fluctuations, the 
commonly used fluctuation quantity 4t flTj is also pre- 
sented in order to compare this measurement to previous 
results l| . The quantity d is related directly to </>t via 

4>T = (CT(T, data) ~ &{T, baseline))**/ < > 

= d<M T > V< N T >, (3.9) 

where Nt represents N tra cks or N c i us . The quantity 4>t 
is related to F T by 

4>T = F T (J(T, baseline) V< N T >■ (3.10) 

The standard deviation of the semi-inclusive spectra can 
be approximated by <T( T ,inci.) ~ <^(T, baseline) V< N T > 
IP^H , where cr(r,incl.) is the standard deviation of the 
semi- inclusive distribution as defined in Eq. 3.3. There- 
fore, 4>t is simply the fraction of non-random fluctua- 
tions in the event-by-event mean px or ex, Ft, scaled by 
(T(r,incl.)- An advantage of Ft over 4>t is that measure- 
ments expressed in Ft can be directly compared without 
further scaling. 

The magnitudes of any non-random fluctuations are 
established by comparing the data distributions to the 
mixed event distributions, which serve as the random 
baseline distributions. For this purpose, the mixed 
event distributions are normalized to minimize the x 2 
value with respect to the data distributions. Figure 
3 and Figure 5 show the M PT and M eT distributions 
for all four centrality classes (data points) with the 
corresponding mixed event M pT and M eT distributions 
overlayed on the data as dotted curves. The broad- 
ening of the distributions for less central collisions are 
due to the reduction in < Ntracks > or < N c i us >. 
Shown in Figure 4 and Figure 6 are the residuals be- 
tween the data and mixed events, defined for each bin i 
as residuak = (M {T ,data), ~ M( T ,mixed)i)/<ri> in units of 
standard deviations, for each centrality class. The shapes 
of the residual distributions are primarily driven by the 
normalization procedure applied to the mixed events. 

For the M PT distributions, the data and mixed event 
distributions are indistinguishable. However, the upper 
M eT edges of the data and mixed event M eT distribu- 
tions show good agreement while the lower M eT edge of 
the data distributions are slightly wider than the mixed 
event distribution. If this low effect were physical, 
it would imply fluctuations with slightly more low eT 
photons since the effect is not seen in the M PT distribu- 
tion for charged particle tracks. However, some of the 
excess fluctuations at low eT can be attributed to the ef- 
fects of cluster merging previously discussed. The mag- 
nitude of this effect has been investigated using a Monte 



Carlo simulation which calculates M eT after reproducing 
the calorimeter cluster separation distribution, the N c i us 
distribution, and the semi-inclusive er distributions from 
the data. The fluctuations in the M eT distribution with 
this effect included in each event are compared to a simu- 
lated mixed event M eT distribution constructed from the 
same generated dataset using the same procedure that is 
applied to the data. In this manner, it is estimated that 
the cluster merging effect contributes an additional Ft 
= 1.5%, 2.1%, 0.9%, and less than 0.01% to the non- 
random fluctuations for the - 5%, - 10%, 10 - 20%, 
and 20 — 30% centrality classes, respectively. The sim- 
ulation confirms that the cluster merging effect signif- 
icantly contributes only to the lower M eT edge of the 
distribution. The remainder of the excess low er fluctu- 
ations is likely due to correlated low energy background. 
GEANT |nj simulations indicate that the primary back- 
ground contribution is produced by low energy electrons 
and muons that scatter off of the pole tips of the central 
arm spectrometer magnet but still pass the cluster se- 
lection cuts. Because of the difficulty in quantifying the 
contribution of background to the excess fluctuations, the 
present M eT data are taken to indicate an upper limit on 
non-statistical fluctuations rather than an indication of 
true non-statistical fluctuations. 

The values of lot, d, Ft, and (j)T for each centrality 
class using the mixed events as the random baseline dis- 
tribution are tabulated in Table III for M PT and Table 
IV for M eT . The errors quoted for these quantities in- 
clude statistical errors and systematic errors due to time- 
dependent variations over the data-taking period. The 
systematic errors are estimated by dividing each dataset 
into nine subsets with each subset containing roughly 
equal numbers of events. For the M PT analysis, the sys- 
tematic errors contribute to 81%, 88%, 76%, and 75% of 
the total error in lj t and 85%, 88%, 80%, and 85% of the 
total error in the variables d, Ft, and 4>t for the — 5%, 
0-10%, 10-20%, and 20-30% centrality classes, respec- 
tively. The corresponding values for the M eT analysis are 
a 67%, 63%, 81%, and 82% contribution to the total er- 
rors in lot, and a 64%, 63%, 81%, and 82% contribution 
to the total errors in d, Ft, and <p t for each centrality 
class. The cluster merging contribution estimates noted 
above are not applied to the values quoted in Table IV. 

IV. DISCUSSION 

Based upon the fluctuation measurements presented 
here, certain fluctuation scenarios in RHIC Au+Au col- 
lisions at -^/SaT/v = 130 GeV are excluded. For this pur- 
pose, we consider two variations of a model that contains 
two classes of events with a difference of effective temper- 
ature, defined as AT = Ti — T±, where Ti is the inverse 
slope parameter of the event class with the higher effec- 
tive temperature, and T\ is the inverse slope parameter of 
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the event class with the lower effective temperature. The 
first variation, Model A, will consider a case where the 
means of the semi-inclusive pt spectra for the two event 
classes are identical, but the standard deviations are dif- 
ferent. The second variation, Model B, will consider a 
case where the means of the semi-inclusive pr spectra are 
different, but the standard deviations are identical. Since 
the semi-inclusive pt distribution is an observed quan- 
tity, the two event classes must be constrained in such 
a way that the mean and standard deviation of the final 
semi-inclusive pr distribution remains constant while the 
effect of the fluctuation manifests itself in the M PT dis- 
tribution. 

The dual event class model is applied to the determi- 
nation of the sensitivity to fluctuations in M PT for the 
— 5% centrality class as follows. Returning to the pre- 
scription outlined in |l5[ ] , the semi- inclusive transverse px 
spectrum can be parameterized by the fv(jpT,P, b) distri- 
bution defined in Eq. 3.1. For both model variations, 
the fraction of events in the event class with the higher 
effective temperature is defined as 



1 _ i-9(l-?)(^) 2 



(^■^ events) class 1 



(-^ events) class 1 H~ events) class 2 



(4.1) 



The pr distribution of the combined sample can then be 
expressed as 

fipr) = qT(p T ,pi,bi) + (1 - q)T(p T ,P2,b 2 ), (4.2) 

where T x = l/bi and T 2 = l/b 2 . 

For Model A, the semi-inclusive px distributions of 
each event class are constrained to have the same mean, 
so we require 



[i=pjb = pi/b 1 = p 2 /b 2 . 



(4.3) 



The variance of the final semi-inclusive pt distribution 
for Model A is constrained by 



1 

P 



q_ 

pi 



P2 



(4.4) 



With these constraints, the choice of a value for q and 
the effective temperature of one event class is sufficient 
to extract the remaining parameters from which sensitiv- 
ity estimates for fluctuations in M PT are obtained. 

For Model B, the semi-inclusive px distributions of 
each event class are allowed to have different means, p\ 
and fj, 2 , so the mean of the total semi- inclusive distribu- 
tion can be expressed as fi = qpi + (1 — q)p 2 - Defining 
a mean shift, A/x, as A// = p 2 — pi, we obtain 



M2 = M + gA/i. 



(4.5) 



Allowing p\ — p 2 and applying the constraint that the 
variances of the two event classes are identical, yields 



Pi 



(4.6) 



With a choice of values for q and Afi, the remaining pa- 
rameters can be calculated, including AT. 

Both variations of the dual event class model arc im- 
plemented in a Monte Carlo simulation in the following 
manner. The number of particles in an event is deter- 
mined by sampling the Ntracks data distribution, approx- 
imated by a Gaussian distribution fit to the data. The 
Pt of each particle in an event is determined individually 
by sampling the appropriate T(j>T,p,b) distribution fit 
to the semi-inclusive pt data distribution, which yields 
p = 0.8 and 6 = 2.46 for - 5% centrality. The p T of 
each particle is restricted to the pt range of the measure- 
ment. With Ntracks and the pt distribution determined, 
the M PT for a given number of events is calculated. The 
generated M PT distribution with q = for either model 
variation is found to be statistically consistent with the 
mixed event M PT distribution. 

The data contain a fraction of background particles 
that did not originate from the collision vertex that effec- 
tively dilute the sensitivity to non-random fluctuations. 
To address this, a fraction of the particles in an event 
are randomly tagged as background particles, whose pt 
distribution is then generated with a separate parame- 
terization prior to calculating M PT for an event. The 
level of background contamination is estimated by pro- 
cessing HIJING Q Au+Au events through a software 
chain that includes a detailed GEANT simulation |nj 
with the complete PHENIX detector geometry included, 
followed by a detailed simulation of the detector elec- 
tronics response [ p"2"[ , whose output is then processed by 
track, cluster, and momentum reconstruction using the 
identical software and input parameters as is used for the 
data analysis. It is estimated that 11% of the tracks and 
26% of the clusters are due to background particles, in- 
dependent of centrality class over the centrality range of 
these measurements. The estimated pt and er distribu- 
tions for the background particles are well parameterized 
by exponential distributions. Again, the majority of the 
er background occurs at low ey, so any correlated back- 
ground would most likely contribute to the lower side of 
the M eT distribution. 

To determine the sensitivity to fluctuations within the 
dual event class model, the fluctuation fraction, g, and 
the value of p\ for Model A and Afj, for Model B are var- 
ied and the M PT distribution is generated at each step. A 
chi-square test is then performed on the generated M PT 
distribution with respect to the mixed event data M PT 
distribution. For a given value of q, the x 2 result in- 
creases as AT increases, which allows a fluctuation exclu- 
sion region to be defined for the single degree of freedom. 
The curves in Figure 7 show the lower exclusion bound- 
aries for the — 5% centrality M PT measurement at the 
95% Confidence Level as a function of q and AT for both 
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variations of the model. If the sensitivity is determined 
based upon the non-mixed data distribution, the lower 
exclusion boundary increases by less than 2 MeV for all 
values of q for either model. Also, for all values of q in 
either model, the estimated background contribution de- 
grades the sensitivity estimates by AT = 3 MeV for both 
models. 

A recent model of event-by-event fluctuations where 
the temperature parameter T = l/b fluctuates with a 
standard deviation ctt on an event-by-event basis [p0| , 
can be simply related to Ft- 



2 Ft 



< T > 2 



p(< n > -1) : 



(4.7) 



where p = 0.8 is the semi-inclusive parameter extracted 
from the present data. For the 0—5% centrality class, the 
present measurement establishes a 95% Confidence Limit 
of 2.6 xl0~ 3 for erf,/ < T > 2 , or 5% for a T / <T>. 
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V. CONCLUSIONS 



The fluctuations in the M PT distributions for all cen- 
trality classes are consistent with the presence of no fluc- 
tuations in excess of the random expectation. The mag- 
nitude of Ft in all cases is positive, which may be due to 
the presence of Hanbury-Brown-Twiss correlations. The 
fluctuations in the M eT distributions do have a small non- 
statistical component, much of which is attributable to 
the effects of merged clusters, the remainder of which are 
taken to indicate an upper limit on non-statistical fluctu- 
ations in transverse energy. By defining a dual event class 
model, limits are set on the amount of M PT fluctuations 
that can be present in the angular aperture \rj\ < 0.35 
and A0 = 58.5° in ^s^n — 130 GeV Au+Au collisions. 
During the RHIC run of 2001, PHENIX has taken data 
for ^/snn = 200 GeV Au+Au collisions with about a 
factor of four increase in azimuthal angular acceptance 
for both the M PT and M eT analyses, which will allow 
the measurements to be extended toward more periph- 
eral collisions. 
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FIG. 1. The Ntracks distribution for the 0— 10% centrality 
class (data points) compared to the N mix distribution from 
the mixed event sample (curve). 
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FIG. 3. The M PT distributions for all centrality classes. 
The curves are the random baseline mixed event distributions. 
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FIG. 2. The M PT distribution for the - 5% centrality 
class. The curve is the result of a 1 distribution calculation 
with parameters taken from the semi-inclusive pr spectra. 



FIG. 4. The residual distribution between the data and 
mixed event M PT spectra as a function of M PT in units 
of standard deviations for all centrality classes. The total 
X 2 and the number of degrees of freedom for the — 5%, 
0-10%, 10-20%, and 20-30% centrality classes are 89.0/37, 
155.7/40, 163.3/47, and 218.4/61, respectively. 
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FIG. 5. The M eT distributions for all centrality classes. 
The curves are the random baseline mixed event distribu- 
tions. The source of differences in the data and mixed event 
distributions are addressed in the text. 
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FIG. 7. The PHENIX sensitivity to non-random fluctua- 
tions in the two variations of the dual event class model that 
are excluded at the 95% confidence level by the M PT analysis 
in the — 5% centrality class. The fraction of events, q, in the 
class of events with the lower inverse slope parameter (event 
class 1) is plotted on the horizontal axis while the difference in 
inverse slope parameter between event class 1 and event class 
2, AT, is plotted on the vertical axis. The curves represent 
the lower boundaries of the excluded regions. 
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FIG. 6. The residual distribution between the data and 
mixed event M eT spectra as a function of M eT in units 
of standard deviations for all centrality classes. The total 
\ 2 and the number of degrees of freedom for the — 5%, 
0-10%, 10-20%, and 20-30% centrality classes are 310.0/32, 
896.4/36, 678.7/47, and 553.9/53, respectively. A large frac- 
tion of the residual contributions are due to the effects of 
cluster merging. 
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TABLE I. Statistics pertaining to the M PT analysis. The 
values of < M PT > are quoted for 0.2 < pr < 1.5 GeV/c and 



are not corrected for efficiency or acceptance. 


Centrality 


0-5% 


- 10% 


10 - 20% 


20 - 30% 


Data 


events 


72692 


149236 


149725 


150365 


^ ^tracks ^ 


59.6 


53.9 


36.6 


25.0 


a N tracks 


10.8 


12.2 


10.2 


7.8 


< M PT > (MeV/c) 


523 


523 


523 


520 


a PT (MeV/c) 


290 


290 


290 


289 


a MvT (MeV/c) 


38.6 


41.1 


49.8 


61.1 


Mixed Events 


< M VT > {MeV/c) 


523 


523 


523 


520 


a MpT (MeV/c) 


37.8 


40.3 


48.8 


60.0 



TABLE III. Fluctuation quantities for the M PT analysis. 



Centrality 



^(T,data)(%) 

d(%) 

F T (%) 

(p PT (MeV/c) 



0-5% - 10% 10 - 20% 20 - 30% 



7.37 ±0.10 7.85 ±0.13 9.52 ±0.14 11.7 ±0.21 

0.14 ±0.15 0.16 ±0.19 0.19 ±0.21 0.21 ± 0.35 

1.9 ±2.1 2.0 ±2.5 2.1 ±2.2 1.8 ± 3.0 

5.65 ±6.02 6.03 ± 7.28 6.11 ±6.63 5.47 ±9.16 



TABLE IV. Fluctuation quantities for the M eT analysis. 



Centrality 



U(T,data)(%) 

d(%) 
F T {%) 
eT (MeV) 



0-5% - 10% 10 - 20% 20 - 30% 



7.32 ±0.07 7.84 ± 0.08 9.58 ±0.17 11.8 ±0.26 

0.30 ±0.09 0.37 ± 0.12 0.38 ± 0.20 0.40 ± 0.32 

4.3 ±1.3 5.0 ±1.6 4.2 ±2.2 3.5 ± 2.8 

11.5 ±3.59 13.6 ±4.23 11.1 ± 5.75 9.28 ± 7.34 



TABLE II. Statistics pertaining to the Me T analysis. The 
values of < M eT > are quoted for 0.225 < e T < 2.0 GeV and 
are not corrected for efficiency or acceptance. 



Centrality 


- 5% 


- 10% 


10 - 20% 


20 - 30% 


Data 


events 


69224 


138882 


140461 


137867 


< N clu3 > 


68.6 


62.1 


41.6 


28.0 




11.6 


13.2 


10.8 


8.3 


< Me T > (MeV) 


466 


462 


448 


439 


ae T (MeV) 


267 


265 


258 


253 


o M , T (MeV) 


34.1 


36.2 


43.0 


51.8 


Mixed Events 


< M eT > (MeV) 


466 


462 


448 


439 


o MllT (MeV) 


32.7 


34.4 


41.3 


50.0 
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